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ABSTRACT 

Weak lensing causes spatially coherent fluctuations in flux of type la supernovae (SNe la). This leasing 
magnihcation allows for weak lensing measurement independent of cosmic shear. It is free of shape 
measurement errors associated with cosmic shear and can therefore be used to diagnose and calibrate 
multiplicative error. Although this lensing magnification is difficult to measure accurately in auto 
correlation, its cross correlation with cosmic shear and galaxy distribution in overlapping area can be 
measured to significantly higher accuracy. Therefore these cross correlations can put useful constraint 
on multiplicative error, and the obtained constraint is free of cosmic variance in weak lensing field. 

We present two methods implementing this idea and estimate their performances. We find that, 
with ^ 1 million SNe la that can be achieved by the proposed D2k survey with the LSST telescope 
(jZhan et al.l 1200811 . multiplicative error of ^ 0.5% for source galaxies at Zs 1 can be detected and 
larger multiplicative error can be corrected to the level of 0.5%. It is therefore a promising approach 
to control the multiplicative to the sub-percent level required for stage IV projects. The combination 
of the two methods even has the potential to diagnose and calibrate galaxy intrinsic alignment, which 
is another major systematic error in cosmic shear cosmology. 

Keywords: Cosmology: the large scale structure: gravitational lensing 


1. INTRODUCTION 

Weak gravitational lensing has great potential of 
probing dark matter, ne utrinos, dark energy and gravity 
at co s mological scales (iRefregied l2003t lAlbrecht et alJ 
20061: iHoekstra fc Jaiii 120081 : iMunshi et al.l l2008t 
Weinberg et al.l 1201311 . All these great applications rely 


on accurate weak lensing measurement. Cosmic shear, 
lensing induced coherent distortion in galaxy shapes, 
can achieve sub-percent statistical error in weak lensing 
measurement. It is therefore a major science driver for 
massive cosmological surveys such as DES, Euclid, HSC, 
LSST and SKA radio survey. However, cosmic shear 
suffers from a variety of systematic errors such as pho¬ 
tometric redshift err ors and galaxy intrinsic alignment 
(for re views, refer to FlSST Scienc e Collaboration et a l. 
( 2009 ): ^_ ILSST D ark Ener^ Science Collaboration 


( 2012 ): iTroxel fc Ishakr ~ (j2014D V. Tremendous efforts 
have been put to understand and correct these system¬ 
atic errors. 

One systematic error which has received intensive 
scrutiny is shear estimation error. It is often con¬ 
veniently decompos ed into a mult i plicat ive error and 
an additive error (iHevmans et al.l I2006D . Stage IV 
projects such as Euclid and LSST put stringent require- 
ment of controlling multiplicative error to ^ 0.2%-0.5% 
(iHuterer et al.l[2?)^ iCropper et al.l [2?)T^ iMassev et al.1 
I201.‘ID . A series of blind community challenges of mas¬ 
sive scale have been car r ied o ut over th e last decade 
STEP l: IHevmans et all (IM)1: STEP2: IMassev et all 
Mil: CREAT08: IBridle et al.l (IMIfil): CREATIO: 
iKitching et al.1 (j2012D : GREATS: iMandelbaum et ^ 
(j2014D L The latest GREATS challenge shows that con¬ 
trol over multiplicative error/bias has been significantly 
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improved. For mock catalogues, various shear estima¬ 
tion methods can control mu ltiplicative error to 1% o r 
even 0.1% when PSF is given (IMandelbaum et'Mll2014l) . 
Further improvement may still be expected by refin¬ 
ing existing shear estimation methods or emerging new 
methods such as the recent l y pro p osed Fourier-space 
method (iZhan g| [Mia IMToLlMn IZhang fc KoiiiMul 
IMIllZhang et al.ll201SD . 

Nevertheless, the performance of shear estimation 
methods depends on many factors such as galaxy size, 
flux (S/N), morphology, selection criteria, weighting 
scheme, a nd the accuracy of PSF interpolation (e.g. 
GREATS: IMandelbaum et all (|2014D L Given unprece¬ 
dented variety of galaxies at z ^ 0—4 of stage IV surveys, 
one must keep caution on whether these shear estimation 
methods can achieve the accuracy estimated from sim¬ 
ulated lensed galaxies. It would then be safer to design 
and apply independent diagnostics of multiplicative bias 
based on real data in a model-independent manner. If 
multiplicative bias is detected by such diagnostics , it can 
then be calibrated consequently. iVallinotto et al.l (I2011D 
proposed to calibrate the multiplicative error against 
lensing magnification in galaxy flux and size, and demon¬ 
strated the potential of such diagnostics. In this paper, 
we propose an alternative method, that is to calibrate 
multiplicative error by lensing magnification of type la 
supernovae (SNe la). 

Supernova (SN) flux is magnified by gravitational 
lensing by a factor /r = 1/[(1 — k)^ — 7^]. Here, fi, k and 
7 i _2 are the lensing magnification/amplification, conver¬ 
gence and shear, respectively. 7^ = 7^ -|- 7I. In the weak 
lensing regime, the measured flux fluctuation of SNe la 
(after standalization) is Sp = fi— l-l-<5™*. Here 5™* is the 
intrinsic flux fluctuation of SNe la. On one hand, this 
lensing magnification contaminates the Hubble diagram 
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and degrades cosmological co nstraints from SNe la 


Friemanlll996l:lHol2ill998l:lDalal et al.ll2003l:ICoorav et alJ 

2006bll. On the other hand, it orovides indeoendent 
measure of weak lensing through the lensing induced 

flux 

fluctuation (IMetca 

1 119991: IHamana & Futamase 

2OOOI: iMenard & Dalai 

20051: iCoorav et al.l l2006al: 

Dodelson & Vallinottol x 

0061: Zentner & Bhattacharva 

200s 

; IBen-Davan & Kalavdzhvanll2014l lAmendola et al. 

2014 

usefi 
allov 
in S 

; lEedeli & Moscardinil l2014ll. and is therefore a 
il source of information. Existing data already 
TS for marginal detection of lensing magnification 
V flux (iKronborg et al.l 120101 iBetoule et al.l 12014 


more z ^ 1 SNe la expected in future surveys, precision 
lensing measurement through SN magnification is very 
promising. 

Lensing measured in this way is free of multiplica¬ 
tive error troubling cosmic shear measurement. The 
measured cosmi c shear 7,— 1.2 can be conveniently 
parametrized as (iHevmans et al.l[2006H . 


7 obs = (1 + +Ci+ 7 “* , ( 1 ) 

with an extra term 7™* arising from the intrinsic galaxy 
shape noise. For the moment we approximate the mea¬ 
sured 7/(1 —k) (reduced shear) as shear 7. The neglected 
complexity will be discussed in 0 7™* in general has a 
dominant component of no spatial correlation, and a spa¬ 
tially correlated component (galaxy intrinsic alignment). 
rrii is the multiplicative error/bias and Ci is the addi- 
tive error. The two mj ca n differ form each other (e.g. 
iMandelbaum et all (|2014D '). For brevity, we will work 
with m defined with respect to cosmic shear E-mode (k). 

The two data sets (cosmic shear and SNe la magnifi¬ 
cation) both measure the same weak lensing, but with 
different prefactors (e.g. 1 -I- m). Therefore combining 
the two data sets we can measure m without assump¬ 
tions on the true lensing signal. Furthermore, if the two 
data sets locate in the same cosmic volume, cosmic vari¬ 
ance of the weak lensing field will be eliminated and will 
not degrade constraint on m. 

The major obstacle in this approach is the low number 
density of SNe la and therefore heavy shot noise. Later in 
this paper we will show that at least half a million SNe la 
are required to diagnose |m| < 1%, the min imum require¬ 
ment for stage IV weak lensing surveys (|Huterer et al.l 
l200b[) . Surveys of a million SNe la with well measured 
light curves are highly ambitious. Nevertheless, surveys 
of such scale can be accessible by telesc opes like the LSST 
telescope and have been proposed (iZhan et al.l 120081 : 
fLSST Dark Energy Science Collaboration! 1201211 . Cos¬ 
mological benefits of such surveys will be many-folds, 
besides the luminosity distance measurement and pecu¬ 
liar velocity measurement. The large scale structure of 
these SNe la allows for measurement of baryon acoustic 
oscillation, which can significantly im prove cosmologica l 
constraints from weak lensing alone (jZhan et al.ll2008h . 
The proposed diagnosis of m is a new bonus of such su¬ 
pernova survey. 

To quantify the capability of diagnosing and calibrat¬ 
ing multiplicative bias with SN magnification, we adopt 
t he baseline survey o f SNe la as the D2k survey proposed 
in iZhan et al.l (|2008f) . This proposed five year survey over 


2000 deg^ by the LSST telescope will result in about 2 
million SNe la with well measured light curves. Among 
them, 0.7 million locate at 0.8 < z < 1.2, one of the 
primary target redshift bins for precision weak lensing 
measurement. For LSST cosmic shear, we assume a to¬ 
tal of 3 billion galaxies over 20000 deg^, with a normal- 
ized redshift distribution n„{z) = exp {—z!z^) !{2z‘i) 
(|Huterer et al.ll20^ IZhan fc KnoxI 120061 1 and z* = 0.4. 
The median redshift is 2.675z* = 1.07. The forecasted 
constraint on multiplicative error calibration is sensitive 
to SN survey parameters, but is very insensitive to cos¬ 
mic shear survey parameters. 

This paper is organized as follows. We discuss two 
implementations (Ml and M2) of diagnosing and cali¬ 
brating multiplicative error combining cosmic shear and 
SN magnification in (J2]and )j3l respectively. We discuss 
and conclude in (jH Some technical details of calculation 
are presented in the appendix. 

2. METHOD ONE 

Method one only uses the two measurements, namely 
cosmic shear and SN magnification, to calibrate multi¬ 
plicative error. For theoretical estimation of the expected 
S/N, it is much more convenient to work in Fourier space 
than in real space. The observable will be (5 _f(£) and 
7°^®(£). ^ is an independent multipole mode. Since we 
are only able to measure lensing magnification of SNe la 
at € < 1000 (Fig. H]), we can treat the lensin g field as 
Gaus sian. The corresponding Fisher matrix is (iTegmarkl 

[TMl 

F«/3=EbTr[C-'(^)C,„(£)C-i(£)C.;3(^)] . (2) 

Here, = d/dX^^p) and A„(,g) is the a(/3)-th parame¬ 
ter to be constrained. C{i) is the covariance matrix for 
the given £ mode, 

rm- f C'^ie)+Ni^ (l+m)C^^(£) \ 

( 3 ) 

C''^, C'*' and are the angular power spectra of /x, 7 
(E-mode) and their cross power spectrum, respectively. 

= 47r/skyCr|’/lVsN is the noise power spectrum in SNe 
la magnification measurement, /sky is the fractional sky 
coverage of overlapping SN survey and cosmic shear sur¬ 
vey. For the D2k survey, /sky = 2000/(4 x ISO^/tt) = 

4.8%. ap = \J is the rms dispersion of flux of 

standardized SNe la. A^sn is the total number of SNe 
la in the given sky area and in the given redshift bin. 
TV''' = 47r/skyCTg/A^7 is the noise power spectrum in cos¬ 
mic shear measurement, is the r.m.s ellipticity, and 
Aly is the total number of galaxies in the given sky area 
and redshift bin for cosmic shear measurement. We will 
take the approximation /x ~ 1 -f 2k. On the other hand, 
the E-mode shear 'Je = k. Therefore = 4C'’' and 

In numerical evaluation of Um throughout the paper, 
we adopt ap = 0.1 and CTs = 0.3. ap quoted here is 
solely the intrinsic scatter cr™*. In reality, it should in¬ 
clude that induced by photo-z error. Photo-z error of 
(Tz increases SN scatter to ap ~ cr™*[l -I- a(tTz/(T“‘)^], 
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Figure 1. The lensing power spectrum (solid line) and corre¬ 
sponding noises (dash lines). The source redshift is 0.8 < Zs < 1.2. 
Shot noises, from large to small, are N^/A in SN magnification, 
in cosmic shear and in the weighted galaxy distribution, 

respectively. Forecast of shot noise targets at th e proposed D2k 
survey of 0.72 million SNe la IjZhan et al.ll20QM) with the LSST 
telescope. 


with a = 2{dlnDL{z)/dz)^ = 2.5 at z = 1. For 
cr“* = 0.1 and cr^ = 0.01(1 + z) (|Zhan et al.ll2008[ i , we 
have ap — 1.1 = 0.11. We then conclude that in¬ 

cluding photo-z error does not significantly change our 
forecast. We focus on redshift bin 0.8 < z < 1.2, in 
which 0.72 X 10® SNe la can be observed by the pro¬ 
posed D2k survey. Fig. [T] plots 7V^/4 and for 
0.8 < z < 1.2. The lensing power spectrum is calcu¬ 
lated using the Limber integral, in which the nonlinear 
matter power spectrum i s calculated using the halofit 
model (jSmith et al.ll2003[l . Due to sparse SN samples, 
each single lensing multipole mode is overwhelmmed by 
shot noise dX 80. However, with 2£A£/sky modes for 
each bin of width A£, we can beat down shot noise by 
a factor Therefore we can still measure the 

lensing power spectrum through SN magnification with 
S/N > 3 at £ - 1000 for M/£ = 0.1. 

In the Fisher matrix analysis, we combine all n 
independent £ modes, which we label as £i (f = 
1, • • • ,n). We take the unknown parameters to be A = 
(Ao, Ai, • • • , A„) = (to, CF'(£i), • • • , with Ao = to. 

The Fisher matrix Fq,^ is calculated and inverted in the 
appendix. We hnd that the error in to is 

2edeUyC'>’\e,) 

CT'(iVM/4-h AT)-k ^ 

Immediately we find a fundamental lower limit for the to 



Figure 2. The Icr constraint on m at 0.8 < z < 1.2 as a function 
of the number of SNe la for the two methods (Ml and M2), and 
(Tm, min? the lower limit of calibrating multiplicative error with SN 
magnification. The vertical dashed line denote A^sn = 0-72 x 10^ 
of the proposed D2k survey with the LSST telescope JZhan et al.l 
l200li . 


calibration, 




^W-l/2 






SN 


( 5 ) 


= 5 X 10"® 




Notice that cr^ = (k^) = J{£)/{2Tr))d£/£. This 
fundamental lower limit corresponds to the limit that 
all other sources of statistical errors vanish and the only 
one left is shot noise in SN magnification. 

This limit can only be achieved under the condition 
and N'^ <C C'*'. The first condition is usually 
satisfied (Fig. [T]) since the galaxy population is much 
denser that the SN population. For example, the num¬ 
ber density of cosmic shear galaxies in a LSST-like sur¬ 
vey at 0.8 < z < 1.2 is 600 times higher than that of 
SNe la even for an ambitious D2k survey, resulting in 
A'’' ~ O.lA^/4. It is for this reason that the constraint 
on the multiplicative error to is limited by SN survey 
conhgurations. 

On the other hand, the second condition A^ <C C'*' 
is violated at £ > 500 (Fig. [T]), reflecting non-negligible 
shot noise per multipole mode in cosmic shear measure¬ 
ment. Therefore in reality we are not to reach the limit 

Fig. [2] shows am as a function of Asn- We find that 
the actual am ^ 3(7^, min (Fig. E]). Nevertheless, the 
calibration accuracy on to can reach am = 8 x 10“® for 
Asn = 0.72 x 10® SNe la expected in the D2k survey. 
This const raint is close to the r equirement of 0.5% on to 
for LSST (jHuterer et al.l[2003l and is therefore encour¬ 
aging. Constraints of to for other redshift bins are shown 





















4 



Figure 3. Relative contribution per logarithmic i bin for 0.8 < 
2 < 1.2 of the proposed D2k SN survey, for the two methods re¬ 
spectively. 

in Table [TJ 

Are there possibilities to further improve constraint on 
TO? Fig. [3] shows (S/N)^, the constraining power per 
logarithmic i bin defined through 



For the method discussed in this section (Ml), the con¬ 
straining power peaks at £ ^ 1000 (Fig. [3]). Contribution 
from t > 1000 is suppressed, since 1000 

(Fig. [ij. Since most contribution comes from relatively 
large scale ^ < 10^, the Gaussian approximation adopted 
through the Fisher matrix estimation is valid. However, 
Fig. [T] shows that the leasing signal peaks I ^ 3000, 
so there are rooms for further improvement. We present 
the second method (M2) to do so. 

3. METHOD TWO 

Method two combines the galaxy distribution available 
in the same survey, together with cosmic shear and SN 
magnification, to improve constraint on multiplicative er¬ 
ror. For a survey like LSST, the highest S/N measure¬ 
ment is for the galaxy clustering (Fig. [T]). Next is cosmic 
shear. SN magnification has the lowest S/N. Therefore 
we can utilize the galaxy-SN magnification cross corre¬ 
lation and galaxy-cosmic shear cross correlation to im¬ 
prove the magnification and cosmic shear measurement. 
Combining the two cross correlations allows better de¬ 
termination of TO. 

The galaxy surface overdensity is 


galaxies in LSST have photo-z information, we can ap¬ 
ply a redshift dependent weighting Wg{z) to improve 
the measurement accuracy of to. We then have two 
measures of cross power spectra. is the measured 
SN magnification-galaxy overdensity cross power spec¬ 
trum. C''*'® is the measured cosmic shear-galaxy over¬ 
density cross power spectrum. The “hat” on top of 
corresponding property (e.g. C) denotes the measured 
quantity with measurement errors. We expect that 
(j~fg = (X _|_ rn)C^^/2. Therefore, we can estimate to 
combining the two measurements and C''*'®, 


Ci^s(^£)l2 

C")'9(£) 


( 8 ) 


The requirement here is that the bin size A£ is sufficiently 
large so error in C^^{£) is small {5C^^ <C C'*'®). The 
expectation value (to) = m + 0{m^). Since |to| <C 1, 
the above estimator is virtually free of systematic bias. 
When taking the ratio, cosmic variance in the galaxy- 
shear correlation cancels that in the galaxy-magnification 
correlation because the two share identical cosmic volume 
and hence identical cosmic variance . Constraint on to 
from a single £ bin, assuming Gaussianity, is 


{C3{£) + m'^){m/A + N^) 

2£A£UyCf^9,^i)/A ^ 

1 /C3 + /Nf^/4: +N'^\ 

2£A£Uyr^ V ' 


Here, = 47r/sky/?V^ is the shot noise in the 

weighted galaxy clustering. = ?Vtot/(kFg) is the 

weighted galaxy number in the given cosmic volume, 
and {Wg) = / ng{z)Wg{z)dz/ f ng(z)dz. Finally we will 
combine all multipole bins to constraint to. 




- 1/2 




( 10 ) 


In Eq. m r = j is the cross correlation 

coefficient between the weighted galaxy distribution and 
leasing. An important step to reduce the calibration er¬ 
ror is to increase r. Due to the large amount of galax¬ 
ies and strong clustering between them, C® at 

£ < 10^. We then have the luxury to weigh these galaxies 
to increase r. Since we have (photometric) redshift in¬ 
formation of galaxies and we have accurate measurement 
of galaxy bias, we can exert a nearly optimal weight¬ 
ing to galaxies such that their mean redshift distribution 
matches that of the leasing kernel. The weighting is 




Wl{z)Ho 
ng{z)bg{z)H{z) ■ 


( 11 ) 


Here, HT is the leasing kernel defined through 


^ f Sgng(z)Wg(z)dz 
® f ng(z)Wg(z)dz 


K = 


5mWL{z) 


dx 

c/Ho • 


( 12 ) 


ng(z) is the mean galaxy redshift distribution and the Wl{z) = {Wl{z, Zs)) is the leasing kernel averaged over 

redshift integral is over the given redshift bin. Since the source galaxy distribution. Wl{z^Zb) is the leasing 
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Table 1 

Sensitivity to multiplicative error for various redshift bins. The 
labels “Ml” and “M2” denote the two calibration methods. The 
number of SN e la is based on the D2k survey proposed in 
IZhan et al.l ll2008lh which plans to cover 2000 deg^ over five years 
with the LSST telescope. The estimation adopts crp =0.1 and 
results for other values of ap should scale by (j^/0.1. 


■2^max] 

[0.6, 0.8] 

[0.8,1.0] 

[1.0,1.2] 

[0.8,1.2] 

[0.6,1.2] 

.^SN 

0.59M 

0.50M 

0.22M 

0.72M 

1.31M 

crm{Ml) 

1.8% 

1.4% 

1.6% 

0.8% 

0.6% 

crm{M2) 

0.9% 

0.7% 

0.8% 

0.5% 

0.4% 


kernel of a single source redshift Zs (e.g. iRefregieJ (|2003[1 1 


f 2^m(l + 

1 - 

x(U 1 
x(zs) 

lo 




After this weighting, we expect r ~ 1. Under this 
limit, the weighted galaxy angular power spectrum C® = 
Cl — c'^. The requirement to achieve r = 1 is negli¬ 
gible stochasticity in galaxy bias, galaxy stochasticity 
will bring r < 1 and therefore degrade constraint on m. 
Fig. [3] shows that most constraining power comes from 
i < 2000, where stochasticity in galaxy distribution is 
insignificant. Therefore we only expect modest degrada¬ 
tion in constraint of m caused by stochasticity. Given 
complexities in modelling galaxy stochasticity, we will 
simply ignore it in this paper and only caution that am 
of method two can be slightly underestimated. 

Under the condition C® iV®'", we can prove that 
Cm = Cm,min- However, Fig. [T] shows that even this 
condition breaks at > 1000, and even worse, C® < TV®'" 
Sit i > 2000. Therefore the actual constraint on m is 
poorer (um > Cm,min)- Numerical results shown in Fig. [5] 
find Um ~ l-5cm,min Over a wide range 10^ < Nqn < 10®. 
Therefore method two can deliver a factor of 2 better 
constraint on m, with respect to method one (3crm,min —t 
l-5cm,min)- For 0.8 < z < 1.2, the constraining power 
peaks at £ ^ 2000, so method two utilizes more lensing 
information than method one whose constraining power 
peaks at £ ^ 1000 (Fig. [3]). This explains why method 
two works better than method one. 

With 0.72 X 10® SNe la at 0.8 < 2 : < 1.2 that the pro¬ 
posed D2k survey can measure, method two can achieve 
Um = 5 X 10~®. This basicall y meets the requirement 
for LSST (|Huterer et al.ll200^ . Apply method two to 
other redshift bins also turns out excellent constraints 
on m (0.4-0.9%, Table[T]). For this purpose, adding a SN 
survey like the proposed D2k to LSST would be highly 
beneficialEl 

4. DISCUSSION 

Method two is superior to method one in many as¬ 
pects. The constraints on m of various redshift bins are 
shown in Table [TJ am of method two is usually a fac¬ 
tor of 1.5-2 smaller than that of method one, showing 
that method two is superior in statistical error. Further¬ 
more, method two is unbiased to the presence of galaxy 

^ Euclid requires \m\ < 0.2% (ICropper et al.ll2013l : IMassev et al.l 
1201.111 . It would require a dedicated SN survey (if any), with twice 
as many SNe la as the D2k survey, to meet the Euclid requirement. 


intrinsic alignment and additive error in cosmic shear 
measurement. Galaxy intrinsic alignment has negligi¬ 
ble contamination to the galaxy-lensing cross correlation 
measurement, since the weighted galaxy distribution has 
vanishing weighting within the source redshift bin. Ad¬ 
ditive error is expected to be uncorrelated with the large 
scale structure and therefore is not expected to bias the 
galaxy-lensing cross correlation measurement. For these 
reasons, the measured m using method two is insensitive 
to neither contamination of galaxy intrinsic alignment 
nor additive error in shear estimation. 

In contrast, method one is susceptible to galaxy intrin¬ 
sic alignment and additive error. The magnification auto 
power spectrum contributes little to constraining multi¬ 
plicative error since it suffers from much larger measure¬ 
ment error (e.g. Fig. [IJ. Therefore method one basically 
interprets the relative difference between the measured 
cosmic shear power spectrum and magnification-cosmic 
shear cross power spectrum as multiplicative error. How¬ 
ever, galaxy intrinsic alignment and additive error con¬ 
taminate the cosmic shear power spectrum in Eq. [3l and 
therefore can cause fake diagnosis of multiplicative error. 

This potential problem in method one can be rendered 
into valuable measurement of galaxy intrinsic align¬ 
ment/additive error, with the help of method two. Basi¬ 
cally, method two determines m by the ratio of galaxy- 
cosmic shear cross correlation and galaxy-magnification 
cross correlation, free of intrinsic alignment/additive er¬ 
ror. Method one measures the combination of m and 
intrinsic alignment/additive error through the ratio of 
magnification-cosmic shear cross correlation and cosmic 
shear auto correlation. With the measured m from 
method two, we can isolate the combined effect of galaxy 
intrinsic alignment and additive error. Therefore, com¬ 
bining method one and method two, in principle one can 
measure galaxy intrinsic alignment/additive error and 
multiplicative error simultaneously. 

So far we have demonstrated the potential of SN mag¬ 
nification in calibrating multiplicative error in cosmic 
shear measurement. There are a number of caveats in 
the proposed calibration. One is the underlying assump¬ 
tion . In reality, fi = 1 + 2k + Zk? + 7 ^ • • •. 

For cosmic shear we actually measure the reduced shear 
g = 7/(1 — k) = 7 + 7 /c -!-•••. These high order terms 
lead to and the induced difference in the two 

properties is of the order ^ 10“®. With the presence 
of these high order terms, measuring m will rely on mod¬ 
elling these terms and therefore rely on cosmology. These 
complexities can be incorporated by simultaneously fit¬ 
ting m and cosmological parameters, which determine 
these high order terms. 

Another potential p roblem is dust e xtinction by in- 
tergalactic gray dust (jGorasanitil l200(Tj) . Such extinc¬ 
tion causes little reddening and therefore can not be effi¬ 
ciently corrected by conventional reddening recipes. The 
induced flux fluctuation is spatially correlated and there- 
fore biases lensing measur ement from SN magnification 
(|Zhang fc Corasanitill200'7l ) . It is a potential problem for 
calibrating multiplicative error with SN magnification. It 
is also a potential problem for len sing measurement bas ed 
on galaxy flux fluctuations (e.g. iSchmidt et al.l (I 2012 D L 
and multi plicative error calibrat ion with galaxy flux fluc¬ 
tuations (IVallinotto et al.l 1201 Ih . Fortunately, in princi¬ 
ple this problem can be alleviated. Gray dust extinc- 
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tion induces galaxy number density fluctuation, whic h 
differs from that induced by lensing (jYang et al.ll201^ . 
Therefore we can infer the dust extinction through galaxy 
clustering and eliminate it in SN magnification. Never¬ 
theless, given large uncertainty in our understanding of 
intergalactic gray dust, it is an important open question 
to pay attention. 

Error in photometry calibration can also potentially 
cause problem. Its calibration error can be spatially cor¬ 
related. It affects both method one and method two. 
For method one, it directly alters SN flux and biases the 
power spectrum measurement by lensing magnification. 
It also affects the number of galaxies in bins of observed 
magnitude. It then indues a correlation between the SN 
flux fluctuation and galaxy number overdensity. Thus it 
can bias the calibrat ion of m using method two0 The 
approach proposed in IVallinotto et all (1201111 by calibrat¬ 
ing multiplicative error with the lensing induced size bias 
is free of both the gray dust extinction problem and the 
photometry calibration problem. So it provides an inde¬ 
pendent and highly complementary approach to diagnose 
and calibrate multiplicative error. 

SNe la are highly complementary to other cosmologi¬ 


cal probes. They not only contribute as standard candles 
and valuable measures of weak lensing. They have al¬ 
ready provid ed robust measurement of peculiar velocity 


at z < 0.05 (IBonvin et al.l 120061 IHaiiffbolle et al.ll2007t 

IWatkins h Feldman 20071: 

IDai et al.l 120111 

and will in 

the future even to z ^ 0.5 

(|Zhanff & Chen 

1200811. With 


millions of SNe la, they can be used as tracers of large 
scale structure to measure baryon acoustic oscillation 
(|Zhan et al.l I2008D . Our work adds a new application 
of SNe la, and a new reason to include survey of mil¬ 
lion SNe la by the LSST telescope or other weak lensing 
facilities. 
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APPENDIX 

DERIVING THE CONSTRAINT ON MULTIPLICATIVE ERROR USING METHOD ONE 
The Fisher matrix (Eq. [2]) can be decomposed into four blocks. 


F = 


M E\ 
G H j 


(Al) 


M = Foo is in fact a single number, 


M = ^iTr[C-i(£,)C,o(£.)C-'(^,)C,o(^*)] = ^ M, . 




E = is a 1 X n matrix, with components 


E, = Fo, = -Tr [C-i(G)C,o(G)C-i(^,)C,*(£*)] 
77 is a n X n diagonal matrix, with diagonal elements 

77,, = ^Tr [C-i(G)C,,(7,)C-i(G)C,,(G)] . 
The inversion of F can be done by block operation, 

(F-')oo = (M-EH-1g)-' . 

Since H is diagonal, we have 

n 

M - EH-^G = ^ [M, - F,77t1g,] . 


2 = 1 


The error in m is then 

O-m = [a 

^ We thank the anonymous referee for pointing out this issue. 


p-llV 2 

ioo 


^ (M, - F,77t1g,) 


.2 = 1 


- 1/2 


(A2) 

(A3) 

(A4) 

(A5) 

(A6) 
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From Eq. |31 we can do the matrix inversion analytically to obtain C We then plug the expression of C, C ^ and 
into the above equations. Finally we obtain 

_ \ ^ _ 

Ct'(( 1 + m)2iVM/4 + iV7) + iV77VA‘/4 j ~ I ^ C'r(Ni^/4 + N't) +N'fNi^/A 




The sum over independent modes (X]i) can be replaced by the integral in the continuum limit. Finally we obtain 

2£d€/,kyC'7.2(£) 7 


C'7(£)(iV'"/4 + N'T) + N'rNf^/4 


(A8) 


This is the most important result for method one, and is used in numerical evaluations shown in Fig. [2 [3] and Table 

[II 
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